Contact mechanics with adhesion: Interfacial separation and contact area 
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We study the adhesive contact between elastic solids with randomly rough, self affine fractal sur- 
faces. We present molecular dynamics (MD) simulation results for the interfacial stress distribution 
and the wall-wall separation. We compare the MD results for the relative contact area and the 
average interfacial separation, with the prediction of the contact mechanics theory of Persson. We 
find good agreement between theory and the simulation results. We apply the theory to the system 
studied by Benz et al.pj involving polymer in contact with polymer, but in this case the adhesion 
gives only a small modification of the interfacial separation as a function of the squeezing pressure. 



With the rapid development of micro/nano electrome- 
chanical devices in last decade, surface forces play a more 
and more important role in modern technology. This is 
due to the increase of the ratio between the number of 
atoms on the surface and that in the volume. When we 
bring two surfaces together, attractive (and repulsive) 
forces act between them, and a non-zero force is often 
required to separate two solid bodies placed in intimate 
contact 0, a phenomenon referred to as adhesion. 

Adhesion manifests itself in many ways. Thus adhe- 
sion on one hand makes it possible for a Gecko to walk 
on the ceilings or run on a vertical wall [J HI- On the 
other hand, adhesion can lead to the failure of micro or 
nano devices, e.g. micro-sized cantilever beams |f|. Thus, 
if it is too long or too thin, the free energy minimum 
state corresponds to the cantilever beam partly bound 
to the substrate, which leads to the failure of the de- 
vice. However, if the surface roughness is increased, the 
non-bonded cantilever state may be stabilized due to the 
decrease of cantilever-substrate binding energy. 

In reality most surfaces are not atomically flat. Even 
if a surface appears flat at low magnification, when we 
study the surface at higher magnification we usually 
observe surface roughness on small length-scale. Simi- 
larly, when two solids with nominally smooth surfaces are 
brought into contact, generally they do not make contact 
everywhere, but at high enough magnification one usu- 
ally observes many non-contact regions. The study of 
the interfacial separation is essential for describing, e.g., 
sealing [7], capillary adhesion Q or optical interference. 

Contact mechanics between solid surfaces is the basis 
for understanding many tribology processes [ilfiol 11, 12| 
such as friction, adhesion, wear and sealing. The two 
most important properties in contact mechanics are the 
area of real contact and the interfacial separation be- 
tween the solid surfaces. For non-adhesive contact and 
small squeezing pressure, the average interfacial separa- 
tion depends logarithmically on the squeezing pressure [3, 
13[, and the (projected) contact area depends linearly on 
the squeezing pressure [3, El, EH- For adhesive contact, 
however, no numerical results have been presented in the 
literatures to test the contact mechanics theory with ad- 



hesion. In this letter, we study the relation between the 
average interfacial separation and the squeezing pressure 
when adhesion is included. We compare the result of MD- 
simulations with a recently developed contact mechanics 

theory[IiE!,E3,El- We find S ood 

agreement with the 
theory, which represents the first test of the theory when 
the adhesive interaction is included in the analysis. 

We review the contact mechanics theory of Persson 
briefly. It can be used to calculate the stress distribution 
at the interface, the area of real contact and the average 
interfacial separation between the solid walls [13. 14 1. In 



this theory, the interface is studied at different magnifi- 
cations Q = L/\ where L is the linear size of the system 
and A the resolution. The wavevectors are defined as 
q = 2tt/X and qi, = 2tt/L so that £ = qjqh- 

Consider an elastic block with a flat surface in adhesive 
contact with a hard substrate with a randomly rough 
surface. Let cr(x, Q denote the (fluctuating) stress at the 
interface between the solids when the system is studied 
at the magnification £. The distribution of interfacial 
stress 



P(a > = (S(cr-a(x ) 0)). 



(1) 



In this definition we do not include the 5(<r)-contribution 
from the non-contact area. 

For perfect (or complete) contact it is easy to show 
that P(a, C) satisfiesQ 



where 
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/(C) = ^E* 2 q L q 3 C(q). 
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Here E* = E/(l — v 2 ) is the effective elastic modulus. 
The surface roughness power spectrum 



C(q) 
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where z = /i(x) is the surface height at the point x = 
(x,y) and where (..) stands for ensemble average. The 



basic idea is now to assume that ([2]) holds locally also for 
incomplete contact. 

To solve ([2]) one needs boundary conditions. If we 
assume that, when studying the system at the lowest 
magnification £ = 1 (where no surface roughness can 
be observed, i.e., the surfaces appear perfectly smooth), 
the stress at the interface is constant and equal to p — 
-Fn /Aq, where i*N is the load and Ao the nominal contact 
area, then P(a, 1) = S(a — p). In addition to this "ini- 
tial condition" we need two boundary conditions along 
the cr-axis. Since there can be no infinitely large stress 
at the interface we require P(cr, () — > as a — > oo. For 
adhesive contact, which interests us here, tensile stress 
occurs at the interface close to the boundary lines of the 
contact regions. In this case we have the boundary con- 
dition P(— <j a , C) = 0, where er a > is the largest tensile 
stress possible. The detachment stress c a (C) depends on 
the magnification and can be related to the effective in- 
terfacial energy (per unit area) 7eff(C) using the theory 
of cracks jl8j 



.(0 



7cff (C)£g \ 1/2 
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where 



7e ff (C)A*(0 = A 7 A*(C 1 )-tf Bl (0> 



where ^4*(C) denotes the total contact area at the mag- 
nification C, which is larger than the projected contact 
area A(Q. U e \{Q is the elastic energy stored at the in- 
terface due to the elastic deformation of the solids on 
length scale shorter than A = L/C, necessary in order to 
bring the solids into adhesive contact (see below). 

From ([2]) it follows that the area of apparent contact 
(projected on the xy-plane) at the magnification £, A(£), 
normalized by the nominal contact area Ao, can be ob- 
tained from 



A(C) 
Ao 



0a(C) 



da P(a,Q 



(3) 



We denote AiQ/Ao = P p (q), where the index p indicates 
that A(Q/Ao depends on the applied squeezing pressure 
p. The area of (apparent) contact at the highest magnifi- 
cation C = Ci gives the real contact area. For the elastic 
energy U e i we usefl9| 



A E* — 
2 



dq q 2 W(q,p)C(q), 



(4) 



QL 



where q^ and q\ are the smallest and the largest surface 
roughness wave vectors, and[l9( 

W{q,p) = P p {q) [0+(l-0)I#(q)] , 

where j3 = 0.4. The equations given above are solved as 
described in Ref. 17 1. 



Let us now consider the (average) interfacial separation 
u as a function of the squeezing pressure p = Fn/Aq. 
Note that as p increases, u decreases and we can consider 
p = p(u) as a function of it. Energy conservation gives [ljjl 



du p(u)Aq = U 



(5) 



where U = U C \ + U a d is the sum of the elastic energy U e \ 
stored at the interface and given by (|4]), and the adhe- 
sional energy C/ ad = — A7v4*(£i). From ([5]) we get 



p(u) 



1 dU(u) 
Aq du 



We can also consider U as a function of p and write 



or 



p(u) 



du = 



1 dU dp 
Aq dp du 



1 dU 
A p dp 



dp 



Integrating from u = (corresponding to p — oo) to u 
(corresponding to the pressure p) gives 



_L_ ^ ^ ldU 

Ao Jr, P d P 



(G) 



which is very convenient for numerical calculations. 

Let us provide some details about the numerical sim- 
ulations. The molecular dynamics system has lateral di- 
mension L x = N x a and L y — N y a, where a is the lattice 
spacing of the block. In order to accurately study con- 
tact mechanics between elastic solids, it is necessary to 
consider that the thickness of the block is (at least) of 
the same order of the lateral size of the longest wave- 
length roughness on the substrate. We have developed a 
multiscale MD approach to study contact mechanics [ioj]. 
Periodic boundary condition has been used in xy plane. 
For the block N x — N y = 400, while the lattice space of 
the substrate b ~ a/0, where <j) = (l + v / 5)/2 is the golden 
mean, in order to avoid the formation of commensurate 
structures at the interface. The mass of the block atoms 
is 197 a.m.u. and the a = 2.6 A. The elastic modulus 
and Poisson ratio of the block are E = 77.2 GPa and 
v = 0.42. For self-affine fractal surfaces, the power spec- 
trum has power-law behavior C{q) ~ q- 2 ( H + 1 ) j where 
the Hurst exponent H is related to the fractal dimension 
Df of the surface via H = 3 — Df, For real surfaces 
this relation holds only for a finite wave vector region 
qo < q < qi- Note that in many cases, there is a roll-off 
wave vector go below which C(g) is approximately con- 
stant. Here g^ = 2n/L, qo = 3gL, gi = 12gz,- go is named 
as roll-off wave-vector. The physical meaning is that by 
choosing go = 3gz, one can obtain a self- average equiva- 
lent to an average over 9 independent samples. In MD 
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FIG. 1: The logarithm (with 10 as basis) of the probability 
distribution of normal stress a (where a is in units of E*) for 
(a) three different adhesion parameters e = 8eo, 4eo and eo, 
referred to as a, b and c respectively, and (b) with 8eo for 
three different pressures p/E* = -0.00532, 0.00832, 0.06265 
denoted by a,b, and c respectively. 



simulations, the substrate is rigid and fractal with frac- 
tal dimension Df = 2.2 and root-mean-square roughness 
/i rms = 10 A. The calculations are carried out under the 
temperature OK. 

The atoms at the interface between block and substrate 
interact with the potential 



V(r) = 4e 



(?)"-(?)' 



where r is the distance between the pair of atoms. The 
parameter e is the binding energy between two atoms 
at separation r = 2 1,/6 ro. In the calculations presented 
below we have used the r$ = 3.28 A and e = eo, 4eo 
and 8eo, where eo = 18.6 meV. By comparing the total 
energy for the surfaces separated with the case where the 
surfaces are in contact at equilibrium (at zero external 
load) we obtain A7 = 0.69 J/m 2 for the case e = eo, and 
4 and 8 times higher interfacial binding energy for the 
other two cases respectively. 

Now, let us discuss how to define contact on the atomic 
scale when adhesion is included. In the absence of ad- 
hesion we have found that the interfacial stress distribu- 
tion gives the most accurate way of deducing the area 
of real contact [20(. When adhesion is included we use a 



FIG. 2: The logarithm (with 10 as basis) of the probability 
distribution of interfacial separations u (where u is in units of 
the root-mean-square roughness amplitude /i ms ) for (a) three 
different adhesion parameters e = 8eo, 4eo and eo, referred to 
as a, b and c respectively, and (b) with 8eo for three different 
pressures p/E* = -0.00532, 0.00832, 0.06265 denoted by a,b, 
and c respectively. 



cut-off length d c to define contact. It is clear from the 
force-distance curve that the only distinctive point is the 
maximum tensile stress 21]. Thus the solids are regarded 
as in contact when the stress increases with separation, 
and separated when the stress decreases with separation. 
In our case, the critical wall-wall distance is d c = 3.68 A. 
When the interfacial separation d < d c , it is defined as 
contact. Otherwise it is non-contact. More discussions 
about d r 



22] 



can be found in Ref. 
The probability distribution of (perpendicular) stress a 
at the interface is shown in Fig. Q] for (a) the adhesion pa- 
rameter e = eo, 4eo and 8eo, and (b) for e = 8eo for three 
different values of the applied stress p (p/E* — —0.00532, 
0.00832 and 0.06265). Note that when e increases, the 
maximum a c of the tensile stress increases roughly pro- 
portional to e. In fact, we expect a c oc e/a, where a 
is of order a bond length. At the same time the maxi- 
mum repulsive pressure increases but weaker than linear. 
The increase in repulsive stress is, of course, due to the 
additionally adhesional load which acts on the block. 

When the applied load or pressure increases, the max- 
imal tensile stress is unchanged, see Fig. DJb). This is 
the expected result since the maximum tensile stress is 




pressure (GPa) 

FIG. 3: Molecular dynamics (square symbols) and theoretical 
(solid lines) results for the relative contact area A/Ao as a 
function of the squeezing pressure (in GPa) In the calculation 
(solid lines) we have used A7 = 0.7, 2.8 and 5.6 J/m 2 for 
the curves indicated by eo, 4eo and 8eo, respectively. The 
theoretical results are obtained from Eq. 
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FIG. 4: Molecular dynamics (square symbols) and theoret- 
ical (solid lines) results for the (natural) logarithm of the 
squeezing pressure p (in units of E*) as a function of the 
average interfacial separation u in units of the root-mean- 
square roughness amplitude h rraB . The theoretical results are 
obtained from Eq. [5] 



associated with breaking the atomic bonds at the edges 
of the contact regions (which can be considered as crack 
tips), and this stress is of course independent of the load. 
However, the maximum repulsive stress increases with 
load, but the effect is quite small considering the large 
change in the load. 

The probability distribution of interfacial separations 
is shown in Fig. [2] In Fig. HJa) we vary the adhesion pa- 
rameter e and in (b) the applied stress p as in Fig. Q] 
When the adhesion increases, the surfaces are pulled 
closer to each other and the distribution of separations 
becomes narrower. Similarly, when the applied pressure 
increases, the separation between the walls decreases. 



The MD simulations (square symbols) are compared 
with the theory (solid lines) for the relative contact area 
A/Aq as a function of the squeezing pressure p (see 
Fig. [3]). In Fig. 0] we show the (natural) logarithm of 
the squeezing pressure p (in units of E* ) as a function of 
the average interfacial separation u in units of the root- 
mean-square roughness amplitude h rms . The theoretical 
results (solid lines) agree very well with the molecular 
dynamics calculations (square symbols). 

Benz et al.[l| have studied the contact mechanics for 
polymer surfaces. They studied the interfacial separation 
as a function of the squeezing pressure, and found an ab- 
solute (average) slope bigger than the one obtained from 
Finite Element calculations [23| for non-adhesive contact 
between the interfacial separation and the log-scale pres- 
sure. For adhesive contact, both theory and MD sim- 
ulations predict a bigger absolute (average) slope than 
the one for non-adhesive contact (see Fig. [4j. However 
for the system studied by Benz et al.[l| we find that the 
adhesional interaction gives only a very small (~ 20%) 
increase in the absolute (average) slope, which cannot 
explain the experimental results by Benz et al. [l[ . In the 
supplementary material}^ we present a critical analysis 
of this point. 

To summarize, we have presented a molecular dynam- 
ics (MD) study of the adhesive contact between elastic 
solids with randomly rough surfaces. We have calculated 
the contact area and the interfacial separation between 
the elastic solids, and compared the results with the pre- 
dictions of a recently developed contact mechanics model, 
which is based on continuum mechanics. Considering the 
uncertainty in how to define the contact area and the in- 
terfacial separation at the atomistic level, and the small 
size of our MD system, the agreement between the theory 
and the MD results is very good. 

We thank U. Tartaglino for many useful discussions. 
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